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Summary 

Multichannel transcranial current stimulation systems offer the possi¬ 
bility of EEG-guided optimized brain stimulation. In this brief technical 
note I explain how it is possible to use transcranial current stimulation tGS 
(which includes tDCS, tACS and tRNS among others^) electric field real¬ 
istic brain models to create a forward “lead-field” matrix and, from that, 
an EEG inverter for cortical mapping. Starting from EEG we show how to 
generate 2D cortical surface dipole fields that could generate the observed 
EEG electrode voltages. The main tool is the reciprocity theorem derived 
by Helmholtz. The application of reciprocity for the generation of a forward 
mapping matrix (lead field matrix as is sometimes known) is well known 
[Rushl969^], but here we will use it in combination with the realistic head 
models of [Miranda2013^] to provide cortical mapping solutions compatible 
with realistic head model tCS optimization. 

I also provide a generalization of the reciprocity theorem [Helmholtz 1853^] 
to the case of multiple electrode contact points and EEG dipole sources, and 
discuss its uses in non-invasive brain stimulation based on EEG. This can 
be used to guide for optimization of transcranial current stimulation with 
multiple channels based on EEG. This, as far as I know, is a novel result. 


^ RufEni G, Wendling F, Merlet I, Molaee-Ardekani B, Mekonnen A, Salvador R, Soria- 
Frisch A, Grau C, Dunne S, Miranda PC, Transcranial cnrrent brain stimulation (tCS): 
models and technologies, IEEE Trans Neural Syst Rehabil Eng. 2013 May;21(3):333-45. 
^Rush and Driscoll, 1969, IEEE Trans. Biom. Eng. 

^Miranda et ah, Nenrolmage 70 (2013) 4858 

Helmholtz, ”Uber einige Gesetz der Vertheilung elektrischer Strdme in korperlichen 
Leitern, mit Anwendnng auf die thierisch- elektrischen Versuche,” Ann. Phys. Chem., ser. 
3, vol. 29, pp. 211-233 and 353-377, 1853. 
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1 Introduction 

The motivation of this note is to develop ways to derive targeted mul¬ 
tichannel transcranial current stimulation (tCS) protocols from EEG data. 
Since both EEG and tGS are mostly cortically relevant, we would like to de¬ 
velop ways to compare the generated electric fields of tGS with EEG dipole 
cortical held sources. 

There are now many groups developing realistic, FEM-based models of 
current propagation in the human head, solving Poisson’s equation. In the 
Starstim software^ we have now integrated such realistic head modeling of 
the generated electric helds by non-invasive multichannel brain stimulation, 
based on the work of [Miranda2013]. 

How can we use such current propagation (electric held) models in to 
infer an EEG dipole forward model? They are certainly related by the 
reciprocity theorem and by the (not-unrelated) fact that are both governed 
by Poisson’s equation. Once we derive a realistic forward model we can 
develop inverse modeling approaches, to go from scalp to source space. 

In this note I concentrate hrst on the problem of cortical mapping, i.e., 
hnding a dipole distribution on the cortical surface that will generate the 
EEG pattern. Although the method described can be used to create volume 
dipole vector held forward operators, given the fact that in our own current 
multichannel systems the number of electrodes is limited to 32, and the issues 
with volume inverse modeling (a very ill-posed problem) it is reasonable to 
work with more constrained models: we hx the locations of the dipoles to 
the cortical surface and their orientation normal to it, rehecting our current 
understanding that EEG sources are mostly cortical and normally aligned 
to the surface (pyramidal neuron populations). Surface mesh points in our 
models on the cortical surface are of the order of 35,000. 

The transformation from voltage space to dipole space is linear. The 
main task is to hnd the linear transformation (the “lead held” matrix). 
What follows is based on the reciprocity theorem. Once we derive a cortically 
mapper, tCS montages can be optimized to target the desired EEG features 
in cortical space. 

Finally, I provide also an extension of this theorem to include the case 
of multiple contact points and current sources. This can be used to guide 
for optimization of transcranial current stimulation with multiple channels 
based on EEG. 


® http://www.neuroelectrics.com 
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2 Reciprocity 

The reciprocity theorem states that there is a relationship between the 
following variables: 

• and jW: is the observed potential difference between two 

points in the scalp due to a lone dipole current source in the brain, 

(in a volume 5V) 

• and is an imposed, stimulation, current between these 

two points in the scalp and the resulting current density in the brain, 

/( 2 ) 

The reciprocity relation is 

or 

Dropping the reference and other indices, we can just write: 

VaIa = -J{x)-E{x)5V (2.1) 

(see the figure below). This is nicely explained by [Plonseyl963®] paper and 
later book. 

With regards to units, if we input uV and we work with 1 mA current 
as a reference, the output will be in nA m. 

The problem we wish to address is: what is the voltage Vab at a point b 
w.r.t. point a due to a lone dipole The answer is: 

^^(1) _ .P') ■ P^')6V 

^ab - 

^ab 

For more dipoles the superposition principle applies, 

-(1) ^(2) 

^^(1) _ Jn ■ En dV 

^ab - (2) 

n ^ab 

where n refers to voxel location ID. 

®Plonsey R., (1963) Reciprocity applied to volume conductors and the EEG. IEEE 
Trans. Biomed. Electron. 10:912 
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Now, in general, the volume corresponding to each dipole will not be the 
same. A more refined equation is 


p-(2) -*(1) 


= 




r(2) 

^ab 


Now, we have computed the solutions to all combinations of electrodes in 
our EEG cap with respect to Cz. Using this expression we can thus compute 
the voltage of any channel location s w.r.t Cz. 

We can now define the algorithm for the forward model. 


1. Put dipole where you want (a 3D vector). 

^ ( 2 ) 

2. Dot product with the electric field E)) produced by a the reciprocal 
montage at a and b using the realistic head model with a unitary 
current (1=1). 

3. Multiply by volume of dipole. 


You now have the potential at site a w.r.t b {Vab = Va — 14). 



n 


In a triangular mesh, each mesh triangle can be associated to its defining 
three mesh points. So if we define the area of a mesh point to be the area 
of all the triangles it touches, and then this area over all mesh points, each 
triangle will show up three times. Thus, to account for dipole area/volume, 
we can properly assign an area to each mesh m point by the formula 

triangles t touched by meshpoint m 

Do this for every single component of the dipole sources you want your 
model to represent (e.g., all those in the cortex) ... 

Algorithmically, we have as many basis functions as electrodes minus 1, 
defining currents from Cz to other points. We need to loop around them, 
and then loop over dipole space (3N) to create the matrix. The output will 
be a 26 X 3N matrix, the famous “lead” or forward-mapping K matrix. 

An option for use here is to focus on the cortical surface and only on 
normally oriented dipoles. This will make the problem much more amenable. 
Interesting, only the normal E fields generated by stimulation are relevant 
according to first order models of neuron electric field interaction. 
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2.1 Inversion with EEG data as constraints 


Once we have the forward mapping K matrix, we need to smooth the prob¬ 
lem for inversion. An approach I’ve used in the past is to minimise the “cur¬ 
vature” of the solution subject to the data constraints. See [Ruffini2002^]. 

Here we do the same, with a different implementation of the Laplacian 
operator which can be used on a mesh. 

A simpler approach we can employ here is to use the (over the surface) 
distance matrix across cortical surface points. The matrix to use for regu- 
larisation would be 

The equation RJ = 0 says the “the value of J at a point is equal to the 
average of points nearby weighted by distance”. We can start with q = 2, 

d{x, x')~‘^ 

T.x^x'd{x,x')-^ 


f{d{x, x' 

'Lx^x' fid{^ 

E.g., we can use 

f{d{x,x')) = exp[-d{x,x')/do] 



i ^ ikj j - \J \ du ^ do 


More generally, we can write 


R(x, x!) = 6{x, x') — 


We can follow the work in [Ruffini2002]. We aim to minimise ” curvature” 
here defined by R subject to the constraints that the EEG data impose 
{V = KJ), which are added using a family of Lagrange multipliers: 

X = ^J^R^RJ + • (R - KJ) 

and the solution to this problem is given by 

J = [R'^R + K^K - K^{KK^)-^KR^R] K^V 

as explained in the aforementioned technical note on arxiv. 

^ “Spherical Harmonics Interpolation, Computation of Laplacians and Gauge Theory” 
at http://arxiv.org/abs/physics/0206007 
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2.2 Inversion with standard spatial regularisation 

A related approach is to start from a functional in which both data and 
curvature constraints appear treated equally. Mathematically this means A 
is input by hand and not solved for. The solution is 

J = [K^K + aR^R]-^K^V 

where the parameter a is input by hand. 

Another option is Tikhonov regularisation (using I instead of R), 

J = [K^K + aI]-^K^V 

where the parameter a is input by hand. 

Implicit referencing to Cz is used. 


Reciprocity theorem (Hemholtz 1853) 

1) EEG: Dipole generates scalp potential: J(x) ^ Va 

2) tCS; Current generates E field: U ^ E(x) 

VaIa = -J{x)-E{x)5V 



Figure 1 - Top: Reciprocity relation symbols and actors 
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Figure 2 - Sample tomographic solutions with 26 and 38 electrodes. On the 
top, source dipole patch and reconstruction. On the bottom, reconstructed 
EEG from a sample subject (20 electrodes). Red colors represent outgoing 
sources, blue ones ingoing sources. 
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3 Application of reciprocity for EEG-guided multichannel 
current brain stimulation (MtCS) 

Here we discuss how we can derive scalp potentials from EEG to design 
optimal stimulation protocols. 

Using our tomographic approach, the EEG data can be inverted to source 
space to define a target map, and then optimize as discussed [Ruffini2014®]. 
That is, from EEG we derive a cortical (or volume) dipole field map, and 
then we can specify that the generated electric fields by MtCS align or 
anti-align with them with specified strengths and so on. The advantage 
of this approach is that it allows for flexible specification of targets and of 
weight maps. Spontaneous EEG or ERP data can be pre-filtered to bands 
of interest or treated in other ways (IGA, etc). 

Here we discuss a simpler approach (although a less flexible one as well). 

Consider the theorem for two dipole sources measured again at a single 
scalp point a: 

Va, la + Va, la = -J{xi) • E{x 2 ) 6 V + -/(xs) • E{x 2 ) 6 V (3.1) 

but this is just the resulting potential, of course, Vala- Hence, it is clear 
that for arbitrary sources, we can generalize the reciprocity theorem a bit, 

Vala = - j dxJ{x)-E{x) (3.2) 

We can see where this is going ... but let us get there faster. 

Using the reciprocity relations we can use the lead matrix to relate a 
dipole source distribution in the brain to scalp potentials, 

Ua = Rax J{X) 

with implicit summation over repeated indices (vector indices on densities 
dropped for simplicity, also summed). The reference point for currents and 
potentials h has been dropped for notational simplicity. 

Here we recall that K = —Eax, is the lead field matrix associated to 
unitary currents, mapping interior to scalp space (thanks to reciprocity). 
Hence (summing over x implicit), 

Va = -EaxJ{x) 

®RufEni et al., Neuroimage 89 (2014) 216225 
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Let’s multiply (and implicitly sum over repeated indices using the Einstein 
convention) both sides by a vector of injected currents (meaning lab)- 

la ba = ~Ia Eax J{x) 

Now, if we inject currents la, the generated electric fields are 

E{x) = Eax la 


Hence, 


IaVa = -E{x)-J{x) 


or, with the explicit summations written out 


^ 4 K = - y E{x) ■ J{x)dx 


(3.3) 


Finally, we can imagine a future current stimulation system which is 
capable of controlling scalp current density. The corresponding equation 
would be 

[ VadA- Ja = - [ E{x) ■ J{x)dx (3.4) 

Ja Jv 

This rather elegant expression says that if you want generated electric 
fields and EEG sources to be correlated, currents and potentials have to be 
anti-correlated. This gives a simple way to determine optimal stimulation 
currents given scalp potential. Make currents and potentials to be as “par¬ 
allel” as possible, given the constraints. E.g., maximize {laVal subject to 
some constraints on maximal current, current sum, total injected current, 
etc. 


We note that this will not optimize things like the electric field magnitude 
at the sources. Rather the specific component of the electric field parallel 
to the sources. 

A limitation of this model-free optimization is that we will not really 
know what size electric fields we are generating, cannot add a weight map 
to work with weighted correlation, etc. But still, give some constraints on 
currents it provides a recipe to optimize currents to EEG sources. This can 
be especially useful in close-loop applications. 
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3.1 Applications 

The reciprocity theorem above can be used for various applications, includ¬ 
ing 


• Online optimization of MtCS from EEG (model-blind optimization 
based on EEG data) 

• Glosed-loop applications from EEG: monitor EEG, create stimulation 
waveform (e.g., of the form I=c V or other optimizations as discussed 
above) to amplify or reduce EEG 

• Experiments involving “Playing back EEG using MtGS currents (I ? 
V), which may actually make some sense based on the above equation. 

• Theoretical analysis purposes (e.g., EEG/tGS electrode density needed, 
see next section) 


3.2 Spatial resolution of tCS 

We can also use this equation to study the issue of how many electrodes 
(density) we need for stimulation. If we assume that the sources lie in a 
2D surface (cortical mapping), for example, it is easy to argue that the 
map from sources to EEG is 1-1, and that, vice-versa, the map from scalp 
currents to cortical space E-fields is 1-1. Thus, in principle, if we could 
measure the full 2D surface with EEG (continuum EEG) we could invert the 
problem for sources, or if we could stimulate over all the entire scalp surface 
(continuum stimulation), we could reproduce any cortical target. Then, the 
needed number of EEG or stimulation electrodes would be transferred to 
the question of what the spatial scale of cortical sources or target is. 

In other words, the density of EEG or tCS electrodes needed for full 
inversion or target in this simplified model depends on the spatial scale 
(autocorrelation length) of the dipole sources on the cortex. We are ignoring 
here the issue of the finite size of electrodes, which act as a spatial filter. 
That is, we assume infinitely small electrodes. 



